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ABSTRACT 

We perform in-depth dynamical modelling of the luminous and dark matter (DM) content 
of the elliptical galaxy NGC 1407. Our strategy consists of solving the spherical Jeans equa¬ 
tions for three independent dynamical tracers: stars, blue GCs and red GCs in a self-consistent 
manner. We adopt a maximum-likelihood Markov-Chain Monte Carlo fitting technique in the 
attempt to constrain the inner slope of the DM density profile (the cusp/core problem), and the 
stellar initial mass function (IMF) of the galaxy. We find the inner logarithmic slope of the DM 
density profiles to be 7 = 0.6 ±0.4, which is consistent with either a DM cusp (7 = 1) or with 
a DM core (7 = 0 ). Our findings are consistent with a Salpeter IMF, and marginally consistent 
with a Kroupa IMF. We infer tangential orbits for the blue GCs, and radial anisotropy for red 
GCs and stars. The modelling results are consistent with the virial mass-concentration relation 
predicted by ACDM simulations. The virial mass of NGC 1407 is log M v ; r = 13.3 ± 0.2 -Mq, 
whereas the stellar mass is log M* = 11.8 ± 0.1 Mq. The overall uncertainties on the mass 
of NGC 1407 are only 5 per cent at the projected stellar effective radius. We attribute the 
disagreement between our results and previous X-ray results to the gas not being in hydro¬ 
static equilibrium in the central regions of the galaxy. The halo of NGC 1407 is found be DM 
dominated, with a dynamical mass-to-light ratio of M/L = 2601^qq Mq/Lq b . However, 
this value can be larger up to a factor of 3 depending on the assumed prior on the DM scale 
radius. 

Key words: galaxies:star clusters - galaxiesievolution- galaxies: kinematics and dynamics 


1 INTRODUCTION 

The idea that the Universe is dominated by cold, non interacting 
dark matter (DM) and by dark energy (A) has been extensively 
tested over the last two decades. Computer simulations make spe¬ 
cific predictions about the properties of DM haloes at z = 0, some 
of which can be tested with the increasing quality of observational 
data. The emerging picture is that, although the ACDM model can 
generally explain observables on large scal es, it may fail when it 
comes to gala ctic or sub-galactic scales (e.g., IWeinberg et al.l20l3l : 
Kr ou 


The modelling of DM in elliptical galaxies is no¬ 
toriously difficult when com p ared to dwarf galaxies or 


200§; 

Kuzio de Narav et al. 200sl: 

Battaglia et al. 2008!: 

Amorisco & Evans 

2014 IColeetal.1 

2 OI 3 Amorisco et al. 

201? Adams et al 

2014b. Unlike dwarf s 

pheroids, ellipticals are 
ur ignorance about the 

baryon dominated at the very centre. O 


initial mass function (IMF) in ellipticals introduces a degeneracy 
between the DM mass and the stellar mass, which limits our ability 
to test ACDM predictions on galactic scales. Moreover, unlike 
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spiral galaxies whose haloes are traced by HI gas (e.g.. lOh et al .1 
1201 lh . the outer regions of elliptical galaxies are notoriously hard 
to observe because they are optically faint. Nevertheless, they can 
be probed with discre t e tracers, such as plane tary nebulae (e.g., 
Ide Lorenzi et alj 20081: Napolitano et al. 20091). globu lar clusters 
(GCs) (e.g. . lDeason?t alj|2012l : ISchuberth^tayj^ljhor diffuse 
tracers, such as hot X-ray gas (e.g., Humphrey et al.ll2006h . Lastly, 
some mass modelling techniques suffer from the mass-anisotropy 
degeneracy, driven by the fact that the orbital an isotropy of stellar 
systems is very hard to infer from the data l iMamon & Lokasl 
120051) . 

It has been shown that the cumulative effect of model de¬ 
generacies on the model outcome can be alleviated by mod- 
elling of multiple dynamical tracers within the same galaxy (e.g,, 
Walker & Penarrubia 2011: Ischuberth et al.ll20ld: iNewman et al.l 


2013i ;iNapolitano et al .120141) . This approach can drastically reduce 


modelling uncertainties by up to a factor of four. 

In this paper, we model the luminous and DM content of 
the elliptical galaxy NGC 1407 using three independent dynami¬ 
cal tracers: the stars, which probe the innermost effective radius of 
the galaxy (R e ), and GCs, which in turn consist of two indepen¬ 
dent (red and blue) subpopulations, out to 10 R e - The blue and the 
red GCs are tho ught to repr esent different stages of galaxy evolu¬ 
tion (Brodie & Straderu006) and are characterized by distinct kine- 


matic and spatial properties (e.g., 

Puzia et al]2004 

; Lee et alj201Cj; 

Woodlev et al.||20ld: Norris et al. 

2012fclPota et alj 

201?). Our strat- 


egy is to build a self-consistent model with minimal assumptions, 
in the attempt to constrain DM parameters that are directly compa¬ 
rable to ACDM predictions. 

Simultaneous modelling of stars a nd GCs has been car¬ 
ried out_Jn_othe^_eni£tical galaxiesjRomanowsk Yj&Kochanekl 


1 200 ll:l Schuberth et all2010jj2012}; Napolitano et al . 201?). For 

galaxy (M 87, Agnello et alj 2014bl) . it has been shown that this 
method can differentiate between cuspy and cored DM profiles. 
The flattening, or steepening, of the DM inner slope is currently 
very relevant for galaxy evolution theories, because it can be related 
to a nu mber of competing physical processes, s uch as DM con¬ 
traction |Gnedin_et_alJ 20041). barvonic fee dback dCole et aPl201ll : 
lDj_Cmtjx>_etaTj[2013|^jMjirtizzi etjdJ l2013h . or self-interacting DM 
1 Rocha et al]2013l;|Peter et alJl2013T). 

Our choice of NGC 1407 is also motivated by the rich dataset 
collected as part of the SAGES L egacy Unif ying Globulars and 
GalaxieS (SLUGGfQ) survey ^Brodie et ali2014h . and by previous 
claims that NGC 1407 has a very elevate d H-b and mass-to-light 
ratio M/L w 800 Mq/Lq (e.g.. lGouldlll993h . in tension with 
scaling relations for galaxy groups 1 Eke et al] |2006l). On the other 
hand, modelling o f X-ray gas in th is galaxy iHumphrev et al.l2006l : 
IZhang et al] 120071 : ISu et al J |20 1 4|) found values of M/L that are 
lower by a factor of two. 

This paper is structured as follows. The dataset and the observ¬ 
ables needed to solve the Jeans equations are discussed in Section 
[2] The dynamical model and the fitting method are both given in 
Section [3] The modelling results from Section [4] are discussed in 
Section [5] where we compare our findings with literature studies, 
computer simulations and mass estimators. Section [6] summarises 
the results of the paper. 


2 DATA 

In this section we discuss the derivation of the observables needed 
to solve the Jeans equations in Section[3] We are interested in the 
one-dimensional radial distribution of stars and GCs, along with the 
velocity dispersion profile for both tracers. The data are products 
the SLUGGS survey and have been partially discussed in a series 


of papers 1 

Proctor et al. 2009:; Foster et al.l2009; Spitler et al. 2012 

Pota et al. 

2013: Arnold et alj2014l). 


2.1 Target and conventions 

NGC 1407 is a massive E0 gal axy at the centre of the dwarf galaxy 
dominated Eridanus A group dBrough et al .1120061) . It shows mod¬ 
erate rotation along the photometric majo r axis and evidence of 
past AGN activity dGiacintucci et alj|2012l) . In this paper, we de¬ 
rive an absolute magnitude of Mb = —21.84 mag and a total 
luminosity of Lb = 8.53 x 10 10 Lq (see H2.2L The galaxy ap¬ 
pears relaxed in optical imaging, but it shows disturbances both in 
X-ray imaging and in velocity space, as we will discuss throughout 
the paper. The galaxy minor-to-major axis ratio is q = 0.95 and 
it does not show any significant variation with radius out to 280 
arcsec dSpolaor et alj|200? : iLi et all 1201 lh . The adopted po sition 
angle and systemic velocity are PA=55 deg dSpolaor et al.l 1 200?) 
and n sys = 1779 km s^ 1 (from NED). 

The distance to NGC 1407 i s very uncertain. Using the GC 
luminosity function, iForbes et~akl d2006l ) found D = 20.9 Mpc, 
whereas surface b rightness fluctuations return distances between 
25 Mpc to 29 Mpc dTonrv et al.l200ll ; ICantiello et al.l200^) . In this 
paper we will assume D — 28.05 Mpc as in lRusli et al.l d2013ah . 
The distance modulus is therefore (m — M) = 32.24 mag. With 
this distance assumption, 1 arcsec corresponds to 0.136 kpc. 

Throughout the paper, we will refer to the galactocentric 2D 
(projected) radius and 3D (de-projected) radius as R and r, respec¬ 
tively. The galactocentric radius is defined as the circularized ra¬ 
dius R 2 = {qX 2 ) + (Y 2 /q), where A' and Y are aligned with the 
galaxy major axis and minor axis, respectively. All quoted mag- 
nitud es w ere reddening corrected according to the dust map of 
ISchlegel et all(ll998lj . 


2.2 Stellar surface brightness 


We combin ed the H-ba n d sur fac e brightne ss profil e from 
HST/ACS dSpolaor et alj 120081: iRusli et all 12013bh. with 
Subaru /Suprime-Cam data in the (/-band presented in Pota et all 
d2013l) . 


The Subaru (/-band surface brightness profile was derived 
with the XVista software packag^. Bright sources and the nearby 
galaxy NGC 1400 were masked out in order to maximize the signal 
from NGC 1407. The (/-band surface brightness profile was renor¬ 
malized to match the B-band profile because these two filters have 
very similar passbands. The wide-field of Suprime-Cam (34 x 27 
arcmin 2 ) allows us to measure the galaxy surface brightness out to 
440 arcsec from the centre. We convert B magnitudes to solar lumi¬ 
nosities by adopting a solar absolute magnitude of Mq b = 5.48 
mag. The resulting luminosity profile is shown in Figure Q] We can 
see that NGC 1407 has a central stellar core and a bumpy luminos¬ 
ity profile outside 100 arcsec. These features can be fitted via mul- 


1 http://sluggs.ucolick.org 


2 http://astronomy.nmsu.edu/holtz/xvista/index.html 
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o.i i. R [kpc] 10. 



10. R [ kpc l 100. 



Figure 1. Radial surface density profiles of our dynamical tracers. Physical 
scales are given on the top axes. The top panel shows the 5-band stellar 
surface brightness data points (grey), and the best fit Sersic law to the data 
(black line). We do not fit the inner core. The bottom panel shows surface 
density data points for blue (blue circles) and red GCs (red squares) along 
with the Sersic best fits. The grey points are the stellar surface brightness 
data arbitrarily rescaled for comparison purposes. Overall, the three dynam¬ 
ical tracers have very different density distributions. 


tiple Sersid (1963) profiles, as performed by Lietal. (2011) and 

iSSafflS — 

For our current purposes, we perform a single Sersic fit to 
the data points in Figure |T| after masking the innermost 2 arc- 
sec. We find a Sersic index n = 4.67 ± 0.15, an effective ra¬ 
dius R e = 100 ± 3 arcsec = 13.6 ± 0.4 kpc, and an intensity 
I(R = R e ) = 3.5 x 10 s Lq arcsec -2 . From the best fit to 
the data, we infer a B-band absolute magnitude of Mb = —21.84 
mag and a total luminosity of Lb = 8.53 x 10 1o Lq. 

We note that our estimate of R e is a factor of two larger com¬ 
pared to values from the literature (e.g.. ISpolaor et al.ll20080 . but 
it is consistent wit h size-l umi nosity sc aling relations from wide- 
held observations l lKormendv et al.ll2009() . Such large values of R e 
have become increasingly common with the advent wide-held pho¬ 
tometry, which is able to det ect fai nt surface brightness features in 
galaxy haloes (e.g., iDon zelli et al. l201 ill . 


2.3 Globular cluster surface density 

We extracted our GC catalogue from Subaru/Suprime-Cam imag¬ 
ing in g ri hl ters. Details on the data acquisition and reduction are 
given in IPota et alj d201 3h . Briefly, GC candidates were selected 
in colour-colour and colour-magnitude space. The resulting GC 



0.85 0.90 0.95 1.00 1.05 1.10 
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Figure 2. Photometric GC properties as a function of (g — i) colour. From 
the top to the bottom we show the projected effective radius Re, the Sersic 
index n and the background level bg computed for different colour cuts at 
(.9 — Ocut- 


Sample 

Re 

[kpc] 

n 

Ne 

[GCs arc min -2 ] 

bg 

[GCs arcmin -2 ] 

Stars 

13.6 ± 0.4 

4.67 ±0.15 

- 

- 

Blue GCs 

47 ±4 

1.6±0.2 

7±1 

1.4±0.2 

Red GCs 

23 ±1 

1.6±0.2 

20±2 

1.4±0.1 


Table 1. Density profile best fit parameters for stars and GCs. A Sersic 
function was fitted to both stars and GCs. For the latter, a colour cut at 
(fl — *) = 0-98 was adopted to separate blue and red GCs. 


colour distribution is clearly bimodal, with a dividing colour be¬ 
tween blue and red GCs at (g — i) « 0.98 mag. This result is in 
agreement with an FIST-based study of the NGC 1407 GC system 
jForbes et al .1120061) . 

The GC surface density was obtained by binning candidates 
brighter than i < 25.5 mag in circular annuli centered on the 
galaxy, and dividing the resulting number by the area of the cor¬ 
responding annulus. Each annulus contained more than 50 GCs per 
bin. Poissonian errors on the surface density were calculated as the 
ratio \/GCs per bin/Area of the annulus. 

The resulting (blue and red) GC surface density profiles were 
fitted with a modified Sersic function of the form: 


Ngc{R) = N e x exp 



+ bg, (1) 


where N e is the nu merical su rface density at R e , b„ = 2n — 1/3 + 
009876n _1 Jciotti & Bertinfl999l ). and bg is the background con¬ 
tamination level, which was assumed to be homogeneous across the 
image. 

The best fit to the GC surface density profiles of blue and red 
GCs are shown in FigureQ] The best fit values are given in TableQ] 
The red GCs are more centra lly concentrated than the blue GCs, as 
found in other galaxies (e.g., |Forte_et al J2005I ; Ba ssino et alj2006l : 

IStrader et all20lil:lForbes et all2012h. 
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2.4 The effect of the colour cut 

We investigate the effect of the adopted colour cut on the inferred 
parameters ( R s and n in particular) from eq|T] We select the bluest 
and reddest GCs with an initial colour cut at (g — i) = 0.98 mag, 
and then we slide the colour cut toward the blue and red wings of 
the colour distribution, thereby clipping the GCs with intermediate 
colour. We fit eq.[T|for each adopted colour cut and we study how 
the photometric parameters vary with this quantity. 

The results are shown in Figure[2] where the x-axis represents 
the adopted colour cut. We can see that the best fit parameters to eq. 
E] are remarkably constant regardless of the colour cut. The Sersic 
index is stable at n m 1.5 ± 0.2, whereas the effective radii of the 
two subpopulations are different, as found in almost all GC systems 
studied so far. 

We experimented with different bin sizes, finding that 30 to 80 
GCs per bin return results consistent with those in Figure[2] where 
60 GCs per bin was adopted. The scatter introduced by varying the 
bin size across this range is included in the uncertainties in Figure[2] 
In light of these results, we decided to perform our mass modelling 
with our fiducial colour cut set at ( g — *)cut = 0.98 mag. This 
cut was applied both to the spectroscopic and to the photometric 
sample. 

2.5 Stellar velocity dispersion 

We combine stellar kinematic data from two telescopes. We use 
ESO 3.6/EFOSC2 major axis lon g-slit data for th e radial range 
between 0 and 40 arcsec from IProctor et al] (20091) . We also use 
multi slit Keck/DEIMOS data from Foster et al. (2015) for the ra¬ 
dial range between 30 and 110 arcsec. 

The root-mean-square velocity dispersion u rms profile was ob¬ 
tained by folding and averaging the data with respect to the galaxy 
centefl Given the lack of apparent rotation in this galaxy (Foster 
et al. 2015), the v rm s is equivalent to the classic velocity dispersion 
a. We will use the quantity u rms throughout the paper. 

The stellar velocity dispersion profile of NGC 1407 is shown 
in the innermost 100 arcsec of Figure [3] The v Tlns profile declines 
smoothly with radius, and it shows a velocity dispersion bump (a- 
bump) between 5 0 < R < 90 arcsec , which is also detected in 
metallicity space (Pastorello et al.ll2014l ). We will discuss how this 
feature affects our results further in the text. 


2.6 Globular cluster velocity dispersion 

The spectroscopic GC sample was derived from the 9 DEIMOS 
masks discussed in lPota et al.i J2013I ). plus one additional DEIMOS 
mask observed on 2013 September 29. The new observations con¬ 
sisted of 4 exposures of 1800 seconds each. The seeing was 0.73 
arcsec. The instrument configuration, data reduction and da ta anal¬ 
ysis were identical to those described in IPota et all (20131 ). From 
this additional mask, we acquired 23 newly confirmed GCs and 16 
GCs duplicates from previous masks. Overall, our spectroscopic 
sample consists of 379 confirmed GCs. 

From this catalogue, we clip 6 GCs with radial velocities de¬ 
viating more than 3<r from the velocity distribution of the N = 20 

3 Following lNapolitano et alj (2009), we estimate the v Tma of long slit data 
as n 2 ms ~ v ro t/2 + cr 2 because we have data only along the major axis. 
This has only a small impact on the final result, given that v TO t <C <r at all 
radii. 


R [kpc] 

0 20 40 60 80 



Figure 3. Velocity dispersion profile for our dynamical tracers. The v; rms 
profile (defined in l|2.6l and l2.51 for stars, blue GCs and red GCs is shown 
in black, blue and red, respectively. Open red points mark the radial bins 
at which velocity dispersion bumps are detected. The black arrow localizes 
the velocity dispersion bump of the stars, which is also shown by black star 
symbols. Blue and red GCs have very different velocity dispersion profiles 
outside 200 arcsec, whereas red GCs and stars are alike. 

closest neighbours (Merrett et al.l!2003l) . This is performed for the 
blue and red GCs separately because these two subpopulations 
have very different velocity dispersion profiles. We also clip 72 
GCs brighter than uj Cen, the brightest star cluster in the Milky 
Way, which we set as the boundary between GCs and ultra com¬ 
pact dwarfs (UCDs). UCDs were excluded from our analysis be¬ 
cause they represent a tracer population which is kin ema tically 
and spatially distin ct from normal GCs (e.g., IZhang et alj|2015l : 
IStrader et alj |20 1 lh . For w Cen, we assume a magnitude Mi = 
— 11.0 mag (i = 21.2 mag). Our final GC catalogue consists of 
153 blue GCs and 148 red GCs. 

We calculate the root mean square velocity dispersion of the 
two GC subpopulations as: 

tVms = Jj - Vsys ) 2 - (A Vi) 2 . (2) 

where Vi is the radial velocity of the i-th GC and Aty is its 
uncertainty. v aya = 1779 km s _ is the systemic velocity of 
N GC 14 07. Uncertainties were derived with the formulae provided 
bv lDanese et al.l (l980l ). We solve eq.[2]for radial bins with 25-30 
GCs per bin. 

Figure [3] shows that the velocity dispersion of the blue GCs 
increases with radius, whereas that of the red GCs decreases with 
radius. The velocity dispersion of the red GCs shows two bumps 
at R ss 150 arcsec and R « 250 arcsec respectively, as found in 
the stellar kinematic profile, although at a different galactocentric 
radius. The contribution of GC rotation is naturally folded into eq. 
[2] although we find that ro tation is s mall, (v ro t/u) < 0.4 for both 
GC subpopulations (Pota et alll2013l) . 


3 MASS MODELLING 

Our dynamical model is based on the study of the velocity moments 
of stars and GCs around NGC 1407. Apart from spherical symme¬ 
try and dynamical equilibrium, we also assume that all dynamical 
tracers are pressure supported, as suggested by the apparent lack of 
rotation both in the stars and GCs. In this context (i.e., small v TO t), 
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the projected velocity dispersion a is very similar to the root-mean- 
square velocity dispersion u 2 ms = cr 2 + v 2 ot . 


3.1 Method 


We want to relate the observed velocity dispersion profiles 
Vrms(R), and the density profiles I(R) of our dynamical tracers, 
to the total mass in NGC 1407. This can be achieved by solving 
the spherical Jean sl j 1 9 1 5i ) equations, with a solution of the form 
dMamon & Lokasll20051) : 

vLs(R) = J^f~ K (ft £) v(r)M(r)y, (3) 


where I(R) is the 2D (projected) density profile of the tracer, v(r) 
is the 3D (de-projected) density profile and M is the total mass 
of the galaxy within the radius r. The parameter (3 = 1 — oj/cr 2 
is the orbital anisotropy which represents the ratio of the velocity 
dispersion in the tangential direction <r 2 to the velocity dispersion 
in the radial direction er 2 . We assume the anisotropy to be con¬ 
stant with radius. The parameter K is a co mpl icated fu nctio n of 
the anisotropy /3 (equation A16 from lMamon & Lokasll2005l ). We 
obtain u(r) by de-projecting the observed luminosity profile of the 
stars I(R), and the observed number density for the GCs Ngc{R), 
respectively. T he de-projection is perform ed via the numerical ap¬ 
proximation of lPrueniel & Simienl d 1996T) . 

The total mass of the galaxy is: 


M(r) = M*(r) + M d (r), (4) 


where M* and are the mass of the stellar and the DM compo¬ 
nent, respectively. 

The stellar mass is obtained by integrating the de-projected 
stellar luminosity profile i /,: 


M*(r) = 4 - 7 tY* 



(5) 


where Y« = ( M/L ), is the stellar mass-to-light ratio, assumed to 
be constant throughout the galaxy. 

We parametr ize the DM den s ity profile w ith a generalized 
Hernquist profile dHemquistlll990l : [zhadl 1 9961) . also dubbed the 
generalized Navarro-Frenk-White profile (gNFW): 


Pd{r) = ps 



1 + 



7-3 


(6) 


where r s and p 3 are the characteristic DM scale radius and DM 
density, respectively. When r <C r s , the DM pr ofile declines a s 
r 1 . When 7 = 1, the DM has a central cusp dNavarro et alJ 19971) . 
whereas when 7 = 0 the DM has a central DM core. When r>tj 
the DM profile declines as r~ 3 , as f o und f or relaxed DM haloes at 
s = 0 in ACDM simulations dNavarro et alJ 19971) . The integration 
of eq.[ 6 ]gives the cumulative DM mass within a given radius: 


M d (r) 


47 T p s r 3 
u ; 






(7) 


where u = 3 — 7 and 2 Fi is Gauss’ hypergeometric function. 


3.2 Maximum likelihood analysis 

We want to find a set of model parameters able to simultaneously 
reproduce (via eq. [3} the empirical velocity dispersion profiles of 
stars, blue GCs and red GCs. This can be achieved by considering 


that the three dynamical tracers are embedded in the same gravi¬ 
tational potential (stars+DM halo). In our model, the galaxy po¬ 
tential is characterized by the parameters: r 3 , p 3 , 7 (which char¬ 
acterize the DM halo) and by T* (which characterizes the stellar 
mass). On the other hand, the orbital anisotropy of stars, blue GCs 
and red GCs (namely P*,Pb,(3r , respectively) is a physical quan¬ 
tity uniquely associated with each dynamical tracer. 

We adopt a maximum-likelihood approach to find our best 
fit solutions. Assuming a Gaussian line-of-sight velocity distribu¬ 
tion, the log-likelihood of obtaining an empirical v rms for one 
particular dynamical tracer k, given a set of model para meters 
S = {r s ,p s , 7 , Y*,/3fc} is (e.g. JWalker & Penarrubiall201 ll) : 


hi Pk (tt rm s k >S7: ) 



+ ln(27rAu r 2 ma!i ) 

( 8 ) 


where N is the sample size, whereas u rm s and Au rms are the em¬ 
pirical velocity dispersion profile of a given dynamical tracer and 
its uncertainty, respectively. The joint log-likelihood function that 
we want to minimize is therefore: 

3 

logp(u 

rms,k is) = £ log Pk (tt r ms,fc | Rk ) (9) 

k = 1 


with 5 = {r s , P*, ^b, Pr}. 

Our model consists of seven free parameters: 
r 3 , p s , 7 , T*, /3», [3b, /3r- We point out that this dynamical 
model could be arbitrarily parametrized, folding in . for example , 
an orbital anisotropy varying with radius (e.g.. IWolf et alj|2010i ). 
However, given the degenerate nature of the problem, we want to 
keep things simple and avoid to over-parametrized the model. 

To efficientl y expl ore the parameter space, we use the emcee 
python module dForeman-Mackev et alJ 12013l ). which performs 
Markov-Chain Monte Carlo (MCMC) sampling of a parameter 
space, given a likelihood function (eq. E and a prior function. 

For the prior function, we adopt uniform priors over the fol¬ 
lowing ranges: logr s = {l,2},kpc logp 3 = {5, 8 }Mq kpc -3 , 
- log(l —/3) = {-1.5; 1}, 7 = {0, 2} and T* = {4, 11.2}T 0jB . 
The prior function is defined such that the likelihood is forced to 
zero outside the above ranges. 

Our prior on r 3 is based on the values predicted for NFW 
haloes covering a wide range of g alaxy masses: from 1 x 10 11 Mq 
to m 2 x 10 13 Mq l lSchaller et alj2014h . We will explore the effect 
of the chosen prior on r s further in H5.51 

The prior on the stellar mass-to-light ratio includes the values 
predicted by stellar populatio n models for NG C 1407 dZhang et al.1 
120071 : iHumphrev et alj|2003) in the case of alKroupal d200ll) IMF 
(T*, s = 4 - 6.5T 0 b ) and a lSalpeteil d 19551) IMF (T*, s = 
6 — 11.2Yq b ). The large uncertainties on Y*,s come from sys- 
tematics in the stellar population models. 

The prior on 7 was chosen to be posit i ve for consistency with 
similar studies (e. g., iBamabe et al.l |20 1 3 : ICappellari et alj|20l3l : 

1 Adams et al.ll20l4) and because inner DM depression with 7 < 0 
is a regime not yet e xplore d by numerical simulations, with a few 
exceptions (e.g., |Pestril2014l) . 

We stress that the above priors were chosen to accommodate 
the values predicted by computer simulations, when available. This 
is co mmon practice for dynamical modelling of stellar systems 
(e.g.. ICappellari et aill20 1 3l : lAdams et al.ll2014l) . and it forces the 
model parameters to stay in the physical regime. However, varying 
the chosen prior range will indeed affect the final results, although 
not strongly, as we will show in 115.51 

The MCMC approach to data fitting is becoming increas- 
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i n>iI y popular for ma ss modelling of galaxies l lMamon et al . 20 ll : 
lAgnello et al.ll2014hl) . and in astrophysics in general. We found by 
experimentation that MCM C is str ongly pref erred over a classic 
“gridded” approach (e.g JSchuberth et alj|2012h . because it allows 
us to efficiently explore a much wider parameter space with reason¬ 
able computational time. 


3.3 Characteristic DM parameters 


We use our MCMC results to infer the posterior distributions of 
some characteristic DM parameters. We compute the virial ra¬ 
dius r v i r , which is the radius where the mean density of the DM 
halo is A v i r = 101 times the critical density p CT it = 1.37 x 
10 2 Mq kpc~ 3 of the Universe. The virial mass and the concen¬ 
tration parameter are, respectively: 


M v 


4 . 

= r A 


virpcri' vir 


and 


vir 

r. 


( 10 ) 


We also estimate the total mass-to-light ratio computed at the 
virial radius (M v i r /L), the DM fraction at 5 R e /"dm( r = 5 R e ) = 
MpM(r)/M (r), and the logarithmi c radial gradient of the mass- 
to-light ratio ANapolitano et al] 20051) : 


V;T 


R e AT 
T» Ar 


(ID 


where R e is the stellar effective radius and the radial range A r was 
set to Ar = 79 kpc based on the radial extent of the fitted data 
points (see Figure[3]l. 


3.4 Remarks on the data fitting 


The substructures detected in the velocity dispersion profiles of 
stars and red GCs suggest that these features are real and not arte¬ 
facts of the data. Interpreting the meaning of such features is be¬ 
yond the scope of this paper. The u-bump s i n NGC 1 407 could 
be the signat ure of major merg ers igchaueret ai] l2014l) . or mi¬ 
nor mergers Sharma et al .11201 it iKafle et all 20121) . In particular, 
ISchauer et al. d2014 ) showed that the <r-bumps arise soon after the 
major merger involving at least one disc galaxy, and that this fea¬ 
tures are time invariant, meaning that the underlying galaxy has 
reached a state of dynamical equilibrium. More importantly, the a- 
bumps show up only in the stellar and GC components, whereas 
the velocity dispersion profile of the DM halo is smooth after the 

merger. _ 

Unlike the dynamical modelling with a lSchwarzschildl dl979l) 
technique, the second-order Jeans equations used in our work are 
not designed to cope with kinematic substructures. In the case of 
non dynamical equilibrium, the generalized Jeans equations can in 
principl e alleviate the effect of substructures on the mass modelling 
results dFalco et al.ll2013l) . On cluster scales, the impact of kine¬ 
matic substruc tures on the inferred DM halo is negligible inside 
2-3 virial radii dFalco et al]2014l) . Also on galactic scales, the kine¬ 
matic substruc tures have been shown to affect the final results up 
to 20 per cent dYencho et al.ll2006l : iDeason et alj|2012l : IKafle et all 

HoH). 


For our specific case, we will assess the effect of substruc¬ 
tures on our results by manually masking the cr—b umps from the 
D rms radial pr ofile, as suggested by lKafle et al.ld2012l) and recently 
performed by lLane et afl d 201 5l ). In the following, we will focus 
on the results for the masked case only (i.e. we do not fit the open 
points and the stellar symbols in Figure[3]l, and we will discuss the 
unmasked case in H5.5I 


The fit to the stellar data was performed using data points out¬ 
side 2 arcsec because our model ca nnot reproduce the r adial vari¬ 
ation of orbital anisotropy found bv lThomas et alj d2014l) within 2 
arcsec. T his also means that the central supermassive black hole 
dRusli et al.l l2013ah is automatically excluded in our model, al¬ 
though the mass of the supermassive black hole is expected to have 
a minimal effect on the inferred DM halo. 


4 RESULTS 


We run 40 MCMC chains (or walkers) with 5000 steps each, which 
explore the parameter space simultaneously. The final acceptance 
rate of MCMC analysis for all our models is always between 35 and 
40 per cent. We discard the first 15 per cent of the chains to account 
for the “burn-in” period, in which chains move from low likelihood 
regions to high likelihood regions in the parameter space. The mod¬ 
elling solutions and the best fits to the data are shown in Figure [4] 
and Figure[5] respectively. 

Figure [4] shows the ID and 2D posterior distributions for our 
seven free parameters. For each parameter, we draw the contours 
containing 68 per cent (lcr) and 95 per cent (2cr) of the posterior 
distribution. The quoted best fit values and uncertainties in Figure 
[4]are the median of each distribution and the la uncertainties (i.e., 
the 16th and 84th percentiles), respectively. 

Figure [4] shows that some of our model parameters are degen¬ 
erate with each other. These degeneracies arise from the interplay 
between mass and anisotropy, and from the uncertain decomposi¬ 
tion of the mass profile into stars and dark matter. For example, 
there is a degeneracy between the slope 7 and the stellar mass-to- 
light ratio T,, such that large val ues o f T.Amply smaller values of 
7, and vice-versa (see also INewman et al1l2013i) . 

Nonetheless, our joint maximum-likelihood analysis allows 
us to put broad constrains on both T* and 7. The best fit stellar 
mass-to-light ratio in NGC 1407 is more consistent with a Salpeter- 
like IMF, although uncertainties propagate into the Kroupa regime. 
This result agrees withi n la with the value T , }B = 6.6to'gTQ B 
found dynamically by iRusli et al] J2013bh . who assumed a log¬ 
arithmic potential, and_with T*,b = 7.6Yq b found dynami¬ 
cally b^ Samurovic 1 20141) . who assumed an NF W halo. Note that 
both IRusli et af ]2013bh and ISamurovIdl d2014l) assumed a differ¬ 


ent parametrization of the DM halo with respect to our generalized 
NFW halo, as we discuss in 114.11 

The inner slope of the DM halo in NGC 1407 is found to be 
intermediate between a flat core and a cuspy (NFW) profile, and 
we cannot discriminate betwee n these two r egimes . A similar con¬ 
clusion was also reached bv lNewman et alj d2013l) in galaxy clus¬ 
ters. Our knowledge of 7 is therefore limited by the uncertainty 
on the stellar mass-to-light ratio: a Kroupa IMF (Y* iB < 6) im¬ 
plies cuspy DM haloes, whereas a Salpeter IMF implies cored DM 
haloes. Nevertheless, our results rule out, at the 2cr level, a very 
cuspy inner DM slope, a s found in M87 (with slope 7 ~ 1.6) 
by lAgnello et al] d2014bt) . The flattening of the DM density pro¬ 
file may be linked to many compe tin g phy sical processes, such 
as baryonic feedba ck dGovemato et al .1 20121) . self-interacting DM 
iRochaeUtl. 2013), or dynamical friction from baryonic clumps 
l El-Zant e t al.ll20 0ll). AGN feedback , known to have occurred in 
NGC 1407 dGiacmtuccj_eta] J|2012|), can indeed produce central 
DM cores dCole et ah 2011 HdI Cintio et alj|2013l : iMartizzi et al.l 
2013). However, this would raise the question on why the mod¬ 
elling of M87, which also hosts an active AGN, supports a central 
DM cusp, rather than a DM core. 
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Figure 4. One and two dimensional projections of the posterior probability distributions of the seven free parameters of our model. Solid and dashed lines 
are the la and 2 a contours of the distribution. The orbital anisotropy for the stars, blue GCs and red GCs is shown in purple, blue and red, respectively. The 
orange and green lines mark the expected value of Y* for a Kroupa IMF and Salpeter IMF, respectively. We show the linear scale radius r s , although the fit 
was performed with the logarithmic equivalent (logr s ). The median values of each distribution and relative lcr uncertainty are quoted on each panel. Note 
that we report the actual value of p and r s , rather than — log(l — P) and log r s , respectively. 


The orbital anisotropy of the blue GCs is found to be very 
tangen tial (/3b ~ —12), as first suggested by iRomanowskv et al.1 
d2009l) . This finding is in contrast to computer simula tions which 
predict that halo particles re side on radial orbits (e.g. jDekel et alj 
booi ls ommer-LarsenllTOObi ). On the other hand, the results for the 
stars (fit ~ +0.4) and red GCs (fin w +0.7) are indeed sugges¬ 
tive of radial orbits. This result supports the idea that stars and red 
GCs have a similar origin. We recall that the value of fit is the aver¬ 
age orbital anisotropy outside the core radius (we only fit the stellar 
d ata out side 2 arc se c). This result also agrees with the results of 
iThomas et al.i ( l2014lf . who found galaxies with central stellar cores 
to have radial anisotropies outside the core radius. Considering the 
strong deviation from isotropy found in our results, the assumption 
of Gaussianity in our model ( f!3.2l) could in principle bias our fits. 

The posterior distributions for the characteristic DM param¬ 


eters defined in 1 )3.31 are shown in Figure [7] The best fit values 
are given in Table [2] The results from this Table are discussed in 
Section[5] We anticipate that our estimate of M v i r /L is consistent 
with the upper limit of M v i r /L ~ 300 Mq/Lq found bv lSu et ail 
d2014l) (although the reader should see lj5]for a full comparison with 
X-ray results). This result supports the suspicion that NGC 1407 is 
not at the centre of an extremely DM dominated group, but it is very 
sensitive to the assumed prior on r s , as we will discuss in 1 )5.51 

Lastly, we show the cumulative mass profile of NGC 1407 in 
Figure [6] The total mass profile was decomposed into DM mass 
Md(r) and stellar mass M*(r) to show the relative contributions 
of these two components. In Figure [6] one can see that the stellar 
mass dominates within the stellar effective radius R e , and also that 
the modelling uncertainties at this radius are minimized. We will 
further come back to this point in 1 )5.41 
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Model 

7 

log M vir 
[Mq] 

r v ir 
[kpc] 

M vir /L c vir 

[Mq/Lq b ] 

/dm(< 5 R e ) 

V;T 

T* 

[Mq/Lq B ] 

log M* 

[AJq] 

gNFW 

NFW 

Cored 

LOG 

LOG 

0 6~*~ 0 ' 4 
u -°-0.4 

= 1 

= 0 

13.34/°/? 

13.30/g™ 

13.38t°-Ji 

13.72/0-17 

13-51/0;“ 

724 till 

705/™ 

750 + ?33 

9 74 + 1 36 

y ‘^+133 
826+™ 

260/iqq 13. 2/35 

240/®g IO. 3 / 1 / 

289/000 18.9/JI 

/?oq+303 

ooz_226 

385/“ 9 

0-83/0;04 

0-85/°;03 

O.80/0;00 

0.79/0;03 

0.84/0;02 

1 n+0.3 
1 ,u - 0.2 

O-9/o 0 ! 2 

0.8/0;? 

1.0/oj 

7 4+ 1 - 5 
' - 4 -2.0 

6 . 1 //? 

8.8/0;® 

9.4/?/ 

= 6.6 

11 7Q+0-08 

11 . <y_o.i3 

ll-7l/o:io 

11.87/0-04 

11 nn+0.03 
ii.yu_Q 05 

11.74 

prior: 10 < r s < 250 kpc 

gNFW 

0 7 +0 ' 4 

U '' —0.5 

1 3 47 + 0.31 

8 OI/ 14 ® 

OC9+369 q z?+7.0 

OOZ-159 y -°-4.0 

0.82/0;05 

1 o +0 - 4 
1 * u - 0.2 

7 fi + 1,5 
‘ - D - 2.0 

11.80/°;°! 

prior: 10 < r s < 500 kpc 

gNFW 

n s+ 0-4 

U - 8 -0.5 

13.59/033 

880/o« 

466/^? B.1±H 

0.82/0;05 

i o +0 - 4 
1 ,u - 0.2 

7 o+1.4 
' - 8 -2.0 

11-82/^3 


Table 2. Solutions to the Jeans equations for different halo parametrizations. The quoted best fit values are the median of the posterior distributions of each 
parameter (defined in the text), with uncertainties representing the 16-th and 86 -th percentile. Solutions are shown for a generalized NFW halo (gNFW), a 
standard NFW halo ( 7 = 1), a cored halo (7 = 0), and a logarithmic (LOG) halo, respectively. The last two columns give the results for a gNFW model when 
the prior on r s is relaxed to the quoted values. 




Figure 5. Best fit to the velocity dispersion profiles. Top and bottom panels 
show the Dmis profiles of stars (in log scale), and blue and red GCs, respec¬ 
tively. The thick and dashed lines are the lcr envelope of the best fit ?; rms 
profile for a gNFW and LOG model, respectively. The two models can fit 
the data equally well. 


4.1 Alternative models 

Our results cannot distinguish between DM cusps and cores in 
NGC 1407. Nonetheless, it is instructive to examine the outcome 
of our model when an NFW (7 = 1) and a cored halo (7 = 0) are 
assumed. By doing so, the NFW solutions can be compared both 
to ACDM predictions and to previous studies of NGC 1407 which 
adopted this form of parametrization (see Section|5j- Therefore, we 
fix 7 in eq.[6]to either 0 or 1, and solve the Jeans equations for the 
remaining six free parameters. 

The DM parameters inferred from this exercise are given in 
Table [2] and shown in Figure [7] Generally speaking, the peaks of 
the NFW and cored distributions (purple and orange histograms in 



Figure 6. Best fit cumulative mass profile of NGC 1407. Dark grey and 
light grey contours show the lcr and 2 <r uncertainties on the total mass pro¬ 
file (i.e., the gNFW model in Table|2j. Black and purple lines show lcr and 
2cr envelopes of the DM component and of the stellar mass component, re¬ 
spectively. The uncertainties on the mass profile are minimized at the stellar 
effective radius R e re 13 kpc. 


Figure [7]( are consistent with each other, and so are the respective 
best fit parameters. However, Figure[7]also shows that a cored halo 
requires more massive and extended DM haloes. Conversely, the 
posterior distributions of the NFW halo are fairly symmetric with 
smaller uncertainties. 

This difference between cored and cuspy solutions is related 
to the degeneracy between 7 and Y*. Imposing a DM core means 
that more luminous matter can be added in the central regions, ex¬ 
plaining the large values of T*, whereas more DM can be spread 
in the outer regions, explaining the large values of M v ; r and /dm- 
As a consequence, the DM halo is more concentrated (larger c v ir). 
Having a more accurate estimate of Y* would help to break this 
degeneracy. 

We now model the DM halo of NGC 1407 using an alternative 
parametrization of the DM distribution, known as the cored loga- 
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Figure 7. Characteristic DM parameters. Panels show the posterior distri¬ 
butions of the virial radius r v i r , virial mass M v i r , total mass-to-light ra¬ 
tio M v i r /L and concentration parameter c v i r respectively. The black and 
black dashed line corresponds to the solutions for a gNFW halo (i.e., 7 free 
to vary) when kinematic substructures are excluded or included in the fit, 
respectively. Purple and orange colours represent the posterior distributions 
for an NFW (7 = 1 ) and cored (7 = 0) halo, respectively. 


rithmic (LOG) model (iBinney & Tremaine! 19871) : 


M d (r) 


1 vpr 3 
G rl+r 2 


( 12 ) 


where r o and vo are the core radius and the asymptotic circular ve¬ 
locity of the halo, respectively. This model has six free parameters: 

r 0 ,v 0 ,Y*,/3*, 0b, Pr¬ 
omt maximum-likelihood analysis returns ro = 44+kpc 
and vo = 483tg® km s 1 . The estimates of (3 are similar to those 
in Figure [4] The characteristic DM parameters for the LOG model 
are given in Table[2] The results for the LOG model are more simi¬ 
lar to those of the cored gNFW solutions, because both these mod¬ 
els have a central DM core. Flowever, the outer slope of the LOG 
profile (r -2 ) allows for more DM in the outer regions compared to 
a gNFW, which declines instead as r -3 . The fit to the i; rms profiles 
of the three tracers is satisfactory. From Figure [5]it can be seen that 
both a LOG model and a gNFW model can fit the data equally well. 

iRusli et a 1.1 <2013bl) also modelled the DM halo of NGC 1407 
with a LOG potential. They found r o = 10.9 kpc and vo = 
340 km s" 1 , which is only marginally consistent with our results. 
A caveat is that Rusli et al. used different data and a different mod¬ 
elling technique, which may bias the comparison with our results. 
Moreover, they focused on a baryon dominated region of the galaxy 
(< 1 -R e ), meaning that their estimate of T,,b = 6 . 61 q'|Tq b 
should in principle be more accurate than our findings. Therefore, 
we fix Y* in our model to the best fit value from Rusli et al., and 
we study how this affects the inferred DM halo parameters. 

We find r'o = 17 ± 4 kpc and vo = 409+ 3 ® km s _1 , with the 
respective DM parameters listed in Table[2] As expected, removing 
the dependency on T* significantly reduces all modelling uncer¬ 
tainties, and it produces results more consistent with our gNFW 
solutions, but still marginally consistent with those of IRusli et al .1 


d2013bl) . In order for the LOG and gNFW models to be consistent 
with each other, the true value of Y* must be somewhere in the 
range Y* : _b = 5 — 8Yq b . This range is more skewed towards a 
Salpeter IMF according to stellar population models (see 1 )3.2t . 


5 DISCUSSION 

In this section we discuss our results in the light of previous mass 
estimates of NGC 1407 and in the context of ACDM predictions. 
Moreover, we want to address earlier findings that NGC 1407 might 
be at the centre of an extremely dark matter dominated group, and 
we want to compare our results to discrete mass estimators. 


5.1 Comparison with previous studies 


We compare our results with literature dynamical models of 
NGC 1407, with the caveat that these were obtained with differ¬ 
ent datasets and modelling techniques. The latter usually involved 
fixing Y*, 7 and /3, whereas we leave these parameters free to vary 
in our model. 

The mass of NG C 1407 has b een modelled via X-rays 


( Humplue^etaLMTOOG 

20091: D as et al. 120101 : 


20091 : Deason et al.| 


Zhang et al. 


¥ K 

Su et al. 120141). GCs (Romano wsky et al. 


Romanowsk^emd. 


2012 ; IS amiirovid 1201 4l) and stellar 


spectroscopy l Saxton & Ferrerasl 20 id : Rusli et all 1201 3bl : 
iThomas et al.l 120141) . Moreover, the dynamics of the dwarf 
galaxies around NGC 1407 suggest a total mass-to-light range 
M^/L sa 3 0 0 - 2500M^/L^ a n d M v j T = 10 13 ~ 1Q 14 M^ 
l Gould 1 19931 : lOuintana et al] 1 19941 : iTullvl 120051 : iTrentham et al.l 
12004 On the other hand, the dynamics of X-ray, stars, and GCs 
are suggestive of a normal DM content for this galaxy, with 


Mvir /L « 300Mq jL q . 

Figure [ 8 ] shows the comparison between our inferred mass 
profile and previous studies. The total mass profiles are shown 
in terms of circular velocity d 2 = GM/r, where we recall that 
M(r) = M*(r) + M d {r). 

The hydrostatic equilibrium equation applied to the hot gas 
in NGC 1407 produces steeply rising mass profiles, along with 
a “kink” in the inner reg ions. These features have been detected 
in bri ght ellipticals (e.g.. iHumphrev et al.ll200d : iNapolitano et al .1 
|20I4|) . but they are not seen in our circular velocity curve in Figure 
[ 8 ] The extreme behaviour of X-ray results may be linked to the X- 
ray g as not being in hydrostatic equilibrium (e.g.. iHumphrev et al .1 
120131) . For the specific case of NGC 1407, the non-equilibriu m sce - 
nario is supported by disturbances in X -ray maps dSu et al.ll2014l) . 
In fact, the resulting v c (r) profile of lSu et alj u2014l ~) agrees with our 
results, probably because they analyzed a relatively rela xed region 
of the X-ray gas outside 40 kpc. The results of IHumphrev et all 
J2006h are marginally consistent with ours, but this is probably 
driven by their adopted stellar mass-to-light ratio (Y*^ = 2 — 
3.6Yq b ), which is a factor of 2 smaller than our fitted value 
(Y= 7 . 4 + 2 '(7© b )• Despite the disagreement between our 
results and X-ray results, the latter find normal mass-to-light ratios 
of M v i r /L ~ 300 Mq/Lq, i n agreement with our resu lts. 

The GC-based results of iRoimanowskv et alj d2009l ) (their GI, 
GR and GT models) imply a very massive DM halo for NGC 1407 
(M v i r = 3.4 — 26 x 10 13 Mq at a distance of 20.9 Mpc), whereas 
we find M v i r = 2.21+ 3 x IO^Mq. Using a GC dataset very 
similar to ours, ISamurov! 1 d20l4 inferred an NFW halo with 
Mvir = 1 x 10 13 Mq, the lowest among the literature studies of 
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Figure 8. Comparison with literature studies and with ACDM predictions. Left Panel. Circular' velocity profiles, expressed as v c = (GMjr) 1 1 . The plot 
shows our best fit circular velocity curve with lcr (grey) and 2a (light grey) envelopes. The black dashed lines are the lcr results for the stellar mass only. 
X -ray results are labelled and coloured accordingly. Uncertainties in most cases are not shown for clarity. The red dashed line correspond to the X-ray analysis 
of iRomanowskv et alj 120091) . The vertical black line marks the outermost data point used in our dynamical model. Right Panel. Virial mass - concentration 
relation. When an NFW halo is assumed, our results are shown as black contours (lcr and 2cr). Results from prev ious studies, when available , are colour coded 
as on the left panel. T he red star and red circle correspond, respectively, to the GC-based and X-ray analysis of I Romanowsky et alj i2009l) . The black star is 
the GC based result of|s amurovi ]{2o3). The black line and lcr scatter is the expected relation in a ACDM Universe (eq. m- 


NGC 1407, and a stellar mass-to-light ratio T*,s = 7.6Tq b , 
which is in marginal agreement with our NFW solution. 

Our model solution may be characterized by the log-slope of 
the total density profile (stars + dark matter), — 7 tot . Inside 1 R e , 
we find 7tot — 1-3, which is intermediate to the gravitational lens- 
ing results for groups and clusters (M 200 ~ 10 14 and ~ 10 15 Mq, 
respectively: INewman et all 120151 ). Over the range 1-4 R e , the 
slope of 7tot — 1.4 is much shallower than the 7tot ~ 2.3 found 
for fast-rotator early-type galaxies (Cappellari et al. 2015). This 
difference is not simply driven by a trend with stellar mass, as sev¬ 
eral of the fast rotators were of comparable mass to NGC 1407, but 
appe ars to refl ect a disti nction in the halo masses. 

[Deason et alj j2012h used a power-law distribution function 
model for the dynamics of t h e tota l GC system, based on the 
data from IRomanowskv et ail i 20091) . They found v c = 323 ± 
20 km s -1 at 39 kpc, which is lower than all the other results at 
the 3cr level (Figure [8). We suspect this is an effect of the re¬ 
striction of those models to mass density slopes of isothermal and 
steeper (7tot ^ 2). For reference, we estimate a DM fraction 
of /dm = 0.83 ± 0.04 at 5 R e , whereas Deason et al. found 
/dm ~ 0.6 — 0.85 depending on the adopted IMF. 

The trend emerging from the comparison with literature stud¬ 
ies of NGC 1407 is that modelling results can be strongly biased 
depending on a number of factors, such as the type of dynamical 
tracer used, the modelling technique adopted and the assumed (or 
inferred) stellar mass-to-light ratio. 


5.2 Comparison with simulations 

We now discuss our results in the context of a ACDM cosmology, 
and compare our findings with properties of relaxed DM haloes at 


2 = 0. Th e ACDM model predicts a well defined M v j r — c v i r rela- 
tion (e.g jBullock et al.ll200ll : IWechsler et a l. 2002; Schalleretal 
2044). Adopting a Planck cosm ology JPlanck Collaboration et al 
2014 h lDutton & Macciol d2014l) found: 


= 13.7 ( 


My 


\ lO n M 0 


(13) 


with a la scatter of 0.1 1 dex at fixed M v i r . 

In order to compare our results with eq. [13] we rely on our 
NFW solutions for a fair comparison. The right panel of Figure [8] 
shows that our best fit NFW halo is consistent with ACDM predic¬ 
tions within lcr. On the other hand, the steeply rising v c (r) profiles 
from X-ray studies tend to produce large DM c once n tratio ns rel¬ 
ative to ACDM predictions. The results of ISu et aD 120141) are in 
better agreement with our results, as expected. The large DM con¬ 


(jHumphrey et al. 2006j; Ettori et al. 

20 id; 

Oguri & Hamana 2011), 

but not all ( Gastaldello et al. 2007 

), X-ray datasets. However, at 
pancy can be explained by in- 

least in galaxy clusters, this discre 


voking baryon physics dFedelilQO 12l : iRasia et al.ll2013l) . 

For an NFW halo, the inferred logarithmic mass-to-light gra¬ 
dient is V 1 T = l.l/n' 4 . which is consis tent with the theoretical 
M* — V;T relation from lNapolitano et al.ld2005 ). In conclusion, if 
we assume that the DM halo in NGC 1407 is an NFW halo, our re¬ 
sults are fully consistent with ACDM predictions for a wide range 
of T*. 


We cannot put strong constraints on the IMF, but we notice 
that our fitted stellar mass-to-light ratio is always T, > 5T 0jfl 
(at the lcr level) regardless of the DM parametrizations. This re¬ 
sult supports the range of Y* values predicted by a Salpeter IMF, 
although it is still consistent with a Kroupa IMF within lcr. A 
Salpeter IMF in NGC 1407 is supported by recent claims that 
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Author 

r * 


VcO*) 




[kmr 1 ] 

Walker et al. (2009) 

1.0 x 

Re 

2.5 (t^rms )oo 

Wolfet al. (2010J 

1.3 x 

Re 

^ ( v rms )°o 

CaoDellari et al. (2013) 

1.0 x 

Re 

2-5 (», 2 ms )n e 

Amorisco & Evans (2011) 

1.7 x 

Re 

3-4 (ilSe 

Agnello et al. (2014a) 

Rm 


Kv ?m a(R<r) 


Table 3. Mass estimators in terms of circular velocity = GM(r*)/r*. 
r* indicates the location of the pinch radius, with v c being the circular 
velocity at that radius. The constants Rm, K and R a are dependent on the 
Sersic index of the dy namical tracer and were extrapolated from Table 1 of 
lAgnello et al. 1 ' (2014a ). 


the IMF is bottom-heavy in galaxies as massive as NGC 1 407 

dCappellari et alj|2012l : IConrov et al.ll2013l : iTortora et al.ll 20 1 4) . 


5.3 Comparison with mass estimators 

A number of analytic techniques have been proposed to infer 
galaxy masses. Their aim is to bypass time-consuming dynamical 
modelling, and determine the total galaxy mass from two observ¬ 
ables: the projected effective radius of a dynamical tracer and a 
single measurement of its velocity dispersion. This provides the 
galaxy mass at one p articular (3D) radi us, dubbed “sweet sp ot” 
Ichurazov et alj|2010h or “p inch radius” ( lAgnello et al.ll2014ah . It 
was shown ([Wolf et alj 20 ld) that at the pinch radius, r+, the uncer¬ 
tainties from modelling assumptions are minimized. This fact was 
exploited to i nfer the presence of DM cusp s and DM cores in dwarf 
galaxies fe.g. JWalker & Penarrubiall201lh . 

Table[3]lists a representative set mass of estimators (expressed 
in terms of circular velocity), as well as the location of the pinch 
radius proposed by different authors. Note that the velocity disper¬ 
sion v rnls to be used in Table[3]can be either the velocity dispersion 
averaged within a certain radius (o^,), or the velocity dispersion 
measured at one particular radius t) r 2 ms (U). 

We insert the observables from stars, blue GCs and red GCs 
into the equations in Table [3] and derive discrete v c values for each 
tracer. We recall that the effective radii of stars, blue GCs and red 
GCs are 13.6 kpc, 47.1 kpc and 23.2 kpc, respectively. The GC ve¬ 
locity dispersion averaged within a particular radius was calculated 
using eq.[2] For the stars, we compute the luminosity-averaged ve¬ 
locity dispersion within a given radius. The velocity dispersion at 
a given radius was inferred from the interpolated v ims profiles in 
Figure [3] The results are shown in Figure [9] 

Taken as a whole, the different mass estimators sample our cir¬ 
cular velocity curve fairly well, and most estimators are consistent 
with our results within 2 a. However, when taken singly, the mass 
estimators are not consistent with each other and they can return 
total masses which disagree up to 1.3 dex at a fixed radius. 

This disagreement is not unexpected given that some mass es¬ 
timators have classically been developed for dynamical modelling 
of dwarf galaxies, and have rarely been tested o n large ellipti¬ 
cals, w ith some exceptions dCannellari et all |20 1 3l : lAgnello et al .1 
l2014bl) . Potential biases may rise from the significant velocity dis¬ 
persion gradients we find in our data (see Figure [3), in contrast to 
the flattish velocity dispersion profiles invoked by mass estimators 
for dwarf galaxies. However, the fact that mass estimators at large 
radii agree with our results is encouraging. Therefore, mass esti¬ 
mators should be used cautiously in elliptical galaxies with non-flat 
velocity dispersion profiles. 


700 
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> 400 


300 


0 20 40 60 80 100 

r [kpc] 

Figure 9. Comparison of mass estimators. The plot shows our best fit cir¬ 
cular velocity curve with lcr (grey) and 2 a (light grey) envelopes. Differ¬ 
ent symbols correspond to the mass estimators labelled on the top left, and 
defined in Table [3] Purple, red and blue symbols correspond to mass esti¬ 
mators applied to stars, red GCs and blue GCs, respectively. Note that our 
modelling uncertainties are minimized at 1 R e . 
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□ Walker et al. ( 2009) 

" • Wolf et al . ( 2010) 

A Cappellari et al . ( 2013) 
OAmorisco & Evans ( 2011) 
-lArAgnello et al . (2014 a) 



5.4 The pinch radius of NGC 1407 

The cumulative uncertainties of our model results are minimized 
at 13.3 kpc, as shown by the black line in Figure QU This ra¬ 
dius is in remarkable agreement with the projected stellar effec- 
tive radius 1 R e = 13 . 6 kpc, as also found by othe r authors (e.g., 
Cannellari et akl 120131 : IWalker & Penarrubiall201 ll : lAgnello et al.l 
2014al) . We find that at 1 R e the uncertainty is only 5 per cent, 
whereas it is 10 per cent between 0.5i? e and 4 R e , and it is never 
larger than 20 per cent in the radial range covered by the data. 

The mass enclosed within lR e is M(< li? e ) = 4.2^2 x 
10 11 Mq. This result is valid for a large range of the orbital 
anisotropies, inner slopes of the DM profile and stellar-mass-to 
light ratios. The DM fraction is foM(lRe) = 0.37j(Qper 
cent, which is similar to the va lues for high-mass galaxies in 
ATLAS 3D dCappellari et al.ll2013l ). Mass estimators at 1 R e give 
M(< 17? e ) = 4.8 ± 0.7 x 10 11 Mq which is also in agreement 
with our results. This confirms that the projected stellar effective 
radius is the optimal radius to measure galaxy masses, albeit not 
the most useful for constraining DM. 

It is interesting to quantify the effect of performing a joint 
multi-tracer dynamical modelling, relative to the results we would 
have obtained by modelling the dynamical tracers independently 
from each other. To do so, we perform dynamical modelling of 
NGC 1407 by using one dynamical tracer at a time (eq. GQ. The 
results are shown as coloured solid lines in Figure [TO] We perform 
the same exercise by jointly modelling the pairs stars-red GCs and 
stars-blue GCs, respectively. These are shown as dashed lines in 
Figure [TOl 

A clear picture emerges from Figure [TO] The modelling of 
GCs produces uncertainties larger than «20 per cent, whereas the 
modelling of the stars is able to constrain M very accurately at R e . 
Outside this radius, the estimate of M from all three dynamical 
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r [kpc] 

Figure 10. Percentage uncertainty on the inferred NGC 1407 total mass 
as a function of radius. Blue and red GCs are colour coded accordingly. 
Solid lines represent the uncertainty on M when stars (purple) and GCs 
are modelled independently. The solid black line is the result for our joint 
multi-tracer modelling. The blue dashed and red dashed lines are the results 
from the joint modelling of stars-blue GCs and stars-red GCs, respectively. 
The joint modelling of blue GCs and stars is enough to reduce uncertainties 
between 5 to 20 per cent in the range covered by the data. 

tracers is very loose. It is only by jointly modelling stars and GCs 
that the uncertainties are drastically reduced down to always < 20 
per cent in the radial range probed by the data. It is worth noting 
that the joint modelling of stars and blue GCs is enough to achieve 
satisfactory uncertainties, with the red GCs playing only a minor 
role. 


5.5 Systematic uncertainties 

We investigate the effect of our modelling assumptions on the final 
results. The fact that the posterior distributions of some model pa¬ 
rameters (for example r s in Figure [4} tend to overfill the allowed 
parameter space is i ndeed a concern, although very commo n among 
similar studies (e.g jBarnabe et al.l2012l : lKafle et all2014l) . This is 
driven by a set of degeneracies which prevent us from determining 
some parameters with satisfactory accuracy. 

A common practice is to choose uniform priors based on re¬ 
sults from computer simulations, which is also our adopted strat¬ 
egy. This ensures that the model parameters stay in the physical 
regime. However, it is instructive to investigate how the choice of 
the priors affects our results. 

We focus on the scale radius r B only, leaving the analysis of 
the full parameter space as a future exercise. We recall that our 
analysis was carried out with the logarithmic prior 10 < r s < 100 
kpc. We relax the priors on r s to 10 < r s < 250 kpc and 10 < 
r s < 500 kpc respectively. The remaining free parameters were left 
unchanged. The results are listed in Table [2] for the gNFW model 
only. 

One can see that main effect of expanding the prior on r s is to 
boost the uncertainties on the DM parameters towards more mas¬ 
sive and more extended DM haloes. The stellar mass is insensitive 
to our new prior choice because the stellar mass contribution is neg¬ 
ligible at large radii. One should also bear in mind that r 3 cannot be 
constrained if this is larger than the radial extent of our data (« 80 
kpc). For reference, we find r s « 50 kpc. 

We tested the effect of the adopted galaxy distance on the fi¬ 
nal results by modelling the galaxy at a distance of 20.9 Mpc (as in 


iRomanowskv et aPl2009l ). instead of the 28.05 Mpc used through¬ 
out this paper (see i )2.1l l. We find that the distance has little im¬ 
pact (always within 1 a) on the characteristic DM parameters, total 
mass-to-light ratio and DM fraction. On the other hand, a closer 
distance returns a stellar mass-to-light ratio (and therefore a stel¬ 
lar mass) that is a factor of 1.3 smaller, Y* = 5.8j)5'i Yq b , 
and a concentration parameter that is a factor of 1.2 larger c v i r = 
15-8±I;g. 

Our analysis was carried out by masking the kinematic sub¬ 
structures detected for the stars and for the red GCs, for the reasons 
explained in ^3.41 Figure[7]shows by how much some DM param¬ 
eters vary when the kinematic substructures are included in the fit. 
The difference between the masked and unmasked models is small 
(well within the ltr errors) and the impact on the results discussed 
in this paper is negligible. The smaller uncertainties found from un¬ 
masked data is due to the inability of the Jeans equation to fit the 
kinematic substructures as discussed in 1 )3.41 


6 SUMMARY 

We have conducted in-depth dynamical modelling of the elliptical 
galaxy NGC 1407. Our approach consisted of solving, simultane¬ 
ously, the spherical Jeans equations for three independent dynam¬ 
ical tracers in NGC 1407: galaxy stars within one effective radius 
(lR e ), and blue and red GCs out to 10i? e - This technique alleviates 
well-known model degeneracies, lowering the final uncertainties. 

Stellar data were constructed by combining long-slit spec¬ 
troscopy with multi-slit spectroscopy from the Keck/DEIMOS 
multi-object spectrograph. We also analysed 153 blue GCs and 148 
red GCs obtained from DEIMOS. Both stellar and GC data are 
products of the SLUGGS survey. The DM was parametrized with a 
generalized Hernquist profile, with an asymptotic outer slope dic¬ 
tated by ACDM predictions. The DM inner slope 7 was left free to 
vary in the attempt to discern whether the DM in NGC 1407 has a 
central cusp or a core (the cusp/core problem). The stellar mass-to- 
light ratio was also left free to vary within priors imposed by sin¬ 
gle stellar population (SSP) modelling. We used a Markov-Chain 
Monte Carlo method to explore the wide 7-dimensional parameter 
space. 

Although we cannot discriminate between a cored or cuspy 
DM profile in NGC 1407, our results are suggestive of a DM central 
slope shallower than a cosmological Navarro-Frenk-White (NFW). 
This result disfavours a steepening of the DM profile due to adi¬ 
abatic contraction, and it may favour a scenario in which DM is 
evacuated from the galaxy centre via some physical mechanism, 
such as baryonic feedback or self interacting DM. Our knowl¬ 
edge of 7 is limited by the uncertainties on the mass-to-light ratio 
Y*. Our estimate of Y* favours a Salpeter IMF. supporting recent 
claims that the IMF becomes bottom-heavy in more massive galax¬ 
ies. Our results are accurate between 5 to 20 per cent in the radial 
range probed by the data, which is a big improvement compared 
to dynamical models with single dynamical tracers. We find that 
our modelling uncertainties are minimized at the stellar effective 
radius, establishing it to be the best radius for measuring galaxy 
masses. 

We confirm that NGC 1407 is surrounded by a DM 
halo with an inferred total mass-to-light ratio of M v i T /L = 
2602JqqMq/Lq. However, the value of M v i T /L increases up to 
a factor of three if the prior on the DM scale radius is relaxed to 
10 < r s < 500 kpc. Therefore, we cannot determine whether or 
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not NGC 1407 is at the centre of an extremely DM dominated group 
as claimed in previous studies. 

We tested that the disagreement between our results and X- 
rays results is due to the X-ray gas not being in hydrostatic equilib¬ 
rium in the inner regions of this galaxy. When a relatively relaxed 
part of the galaxy is considered, our findings are in better agreement 
with X-ray results. When the DM halo ofNGC 1407 is assumed to 
follow a cosmological NFW profile, our results are consistent with 
the predicted virial mass-concentration relation from ACDM. We 
compared our mass profile of NGC 1407 with a set of discrete mass 
estimators, finding a marginal agreement. However, we argue that 
mass estimators should be used cautiously in galaxies with strong 
velocity dispersion gradients because their results may be strongly 
biased. 
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